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ABSTRACT 

Context. Phase-referencing is a standard calibration technique in radio interferometry, particularly suited for the detection of weak 
sources close to the sensitivity limits of the interferometers. However, effects from a changing atmosphere and inaccuracies in the 
correlator model may affect the phase-referenced images, and lead to wrong estimates of source flux densities and positions. A 
systematic observational study of signal decoherence in phase-referencing and its effects in the image plane has not been performed 
yet. 

Aims. We systematically studied how the signal coherence in Very-Long-Baseline-Interferometry (VLBI) observations is affected by 
a phase-reference calibration at different frequencies and for different calibrator-to-target separations. The results obtained should be 
of interest for a correct interpretation of many phase-referenced observations with VLBI. 

Methods. We observed a set of 13 strong sources (the S5 polar cap sample) at 8.4 and 15 GHz in phase-reference mode with 32 
different calibrator/target combinations spanning angular separations between 1.5 and 20.5 degrees. We obtained phase-referenced 
images and studied how the dynamic range and peak flux-density depend on observing frequency and source separation. 
Results. We obtained dynamic ranges and peak flux densities of the phase-referenced images as a function of frequency and separation 
from the calibrator We compared our results with models and phenomenological equations previously reported. 
Conclusions. The dynamic range of the phase-referenced images is strongly limited by the atmosphere at all frequencies and for 
all source separations. The limiting dynamic range is inversely proportional to the sine of the calibrator-to-target separation. Not 
surpriseingly, we also find that the peak flux densities decrease with source separation, relative to those obtained from the self- 
calibrated images. 

Key words. Techniques: interferometric - Atmospheric effects 
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1. Introduction 

Phase-referencing is a standard calibration technique in radio 
interferometry. It allows the detection of a weak source (target 
source) by using quasi-simultaneous observations of a closeby 
strong source (calibrator) (see e.g., Ros ,2005 1 and references 
therein). This technique also allows the user to recover the 
position of the target source relative to that of the calibrator; 
a position otherwise lost by the use of closure phases in the 
imaging. Basically, phase-referencing consists of estimating the 
antenna-based complex gains with the calibrator fringes, time- 
interpolating these gains to the observations of the target, and 
calibrating the visibilities of the target with the interpolated 
gains. Therefore, it is assumed that for each antenna the optical 
paths of the signals from both sources are similar. However, at- 
mospheric turbulences and/or inaccuracies in the correlator geo- 
metrical model may introduce errors in the estimates of the opti- 
cal paths of the signals and severely affect the phase-referencing. 
These errors can be partially corrected by applying the self- 
calibration algorithm (see, e.g., Readhead & Wilkinson 1 19781 1 
after phase-referencing. However, self-calibration, specially on 
observations of weak sources, may affect the resulting images in 
undesirable ways (see, e.g., Marti- Vidal & Marcaide l2008ll . 

The correlator model includes contributions from the dry 
troposphere, the Earth orientation parameters, gain corrections 



Alexander von Humboldt Fellow 



for the sampling, and feed rotation of the alt-azimuthal mounts 
of the antennae. An imperfect modelling in any of these con- 
tributions and the loss of coherence of the radio waves within 
the time elapsed between consecutive observations of a given 
source, have an impact on the quality of phase-referencing. 
Some authors have stated that the loss of signal coherence in 
phase-referencing is linearly dependent on the separation be- 
tween the calibrator source and the target source (Beasley & 
Conway , 1995) . However, in a recent publication Marti- Vidal et 
al. (1201 Oi l suggest a phenomenological model different from that 
of Beasley & Conway ( 1995 ), based on Monte Carlo simulations 
of atmospheric turbulences. 

To empirically establish this dependence, we need to com- 
pare phase-referenced images to those obtained from the self- 
calibrated visibilities (i.e., the images obtained by applying the 
optimum phase gains; those that maximize the signal coherence 
of the target sources). For this purpose, both calibrator and tar- 
get must be strong enough to generate fringes with a signal-to- 
noise ratio (SNR) high enough to allow for an accurate estimate 
of the phase, delay, and rate of the fringe peaks, thus avoiding 
the bias effects related to self-calibration of weak signals (e.g., 
Mai-ti'- Vidal & Marcaide l2008] l. 

The S5 polar cap sample (a subset of 13 sources from the 
S5 survey, see Kiihr et al. 1981 , Eckart et al. 1986) is an ideal 
set of sources to perform such a study. The spectra of these 
sources are relatively flat at radio wavelenghts and their flux 
densities range from hundreds of mJy to several Jy. Therefore 
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it is possible to study the loss of phase coherence as a func- 
tion of observing frequency and source separation. We observed 
the S5 polar cap sample at 8.4 GHz and 15 GHz. These observa- 
tions are part of a large astrometry programme (Ros et al. 120011 
Perez-Torres et al. .200 4, Marti- Vidal et al. 2008 and Jimenez- 
Monferrer et al. in preparation) and were performed in phase- 
reference mode with many different calibrator/target combina- 
tions. From these observations, we studied the loss of signal co- 
herence in phase-referenced observations by comparing phase- 
referenced and self-calibrated images for all possible calibra- 
tor/target combinations allowed by the observations. 

In the next section we describe our observations and the pro- 
cess of data calibration and reduction. In Sect.|3]we report on the 
results obtained. In Sect. |4] we summarize our conclusions. 

2. Observations and data reduction 

We observed the 13 sources of the S5 polar cap sample with 
the complete Very Long Baseline Array (VLB A) at 8.4 GHz on 
2001 February 3 and at 15 GHz on 2000 June 15. At each epoch, 
8 bands of 8 MHz bandwidth each were recorded, obtaining a 
synthesized bandwidth of 64 MHz in one polarization (LCP). All 
three epochs consisted of 24 hr of observations. The observations 
took place in phase-referenced mode with different subsets of 
three or four radio sources in different duty cycles (see Martf- 
Vidal et al. 120081 for more details on the observing schedule). 
The sources of each subset were observed cyclically for about 
two hours. Each radio source was observed a total of about four 
hours. Data were cross-correlated at the Array Operation Center 
of the National Radio Astronomy Observatory (NRAO). Details 
on the source images at 8.4 GHz and 15 GHz can be found in 
Ros et al. (120011) and Perez-Torres et al. (I2004I I. respectively. The 
results of the high-precision astrometry analysis at 15 GHz can 
be found in Martf-Vidal et al. (2008). Those at 8.4 GHz will be 
published elsewhere (Jimenez-Monferrer et al. in preparation). 

The data calibration and reduction was performed using 
ParselTongue (Kettenis et al. 120061 ). a Python interface to the 
NRAO Astronomical Image Processing System (aips). We gen- 
erated a script in ParselTongue to automatize the calibration and 
imaging of the phase-referenced images between all source pairs 
that were observed within the same duty cycles. A total of 32 
source pairs were obtained, covering a range of source separa- 
tions between 1.5 and 20.5 degrees. We checked the quality of 
the automated calibration and imaging by re-generating man- 
ually phased-referenced images of some pairs of sources and 
comparing them to the automated images. The results obtained 
were compatible within 0.1 cr. Below we summarize the process 
of calibration followed in our analysis. 

- Step 1. A correction of the the Earth Orientation Parameters 
(EOP) as estimated by the United States Naval Observatory 
(USNO) was applied to the data. 

- Step 2. The visibility phases were aligned between the sub- 
bands, through the whole 64 MHz band for all sources and 
times by fringe-fitting the sub-band delays of a selected scan 
with high-quality fringes and applying the estimated delays 
and phases to all visibilities. 

- Step 3. A second fringe fitting, using now the delays deter- 
mined from the whole band, provided new phase corrections 
for all the observations. The visibility amplitude calibration 
was performed with the system temperatures and gain curves 
from each antenna. 

- Step 4. The calibrated data were exported into the pro- 
gramme DiFMAP (Shepherd et al. 11995) to obtain images of 



all sources. The CLEAN algorithm and several iterations of 
phase and gain self-calibration were applied to each source 
until high-quality images (with residuals close to the thermal 
noise) were obtained. 
- Step 5. These images were imported back into aips for a sec- 
ond fringe fitting, now taking into account the contribution 
of the source structures in the estimates of the residual de- 
lays and phases. The amplitude calibration was also refined 
by estimating the amplitude gains based on the source struc- 
tures (one gain solution for all antennae every 10 minutes). 

Once the data were calibrated as described above, we pro- 
ceeded with the analysis. For each pair of sources (A and B), four 
sets of data were generated. On the one hand, the self-solutions 
of A (B) in the time range when this source was observed in the 
same duty cycles as B (A). On the other hand, the set of vis- 
ibilities of A (B) phase-referenced to B (A), obviously for the 
same time range as that of the corresponding self-calibrated data 
sets. Each of these data sets was used to generate an image us- 
ing natural weighting of the visibilities (to optimize the sensi- 
tivity) and the same CLEAN windows that were used to obtain 
the images in Step 4 of the data calibration procedure. For the 
phase-referenced images, the position of the peak flux-density 
was first estimated and the CLEAN windows were accordingly 
shifted before CLEANing. 

3. Results and discussion 

3.1. Dynamic range 

The dynamic ranges of the images obtained as described in the 
previous section were computed as the peak flux-density per unit 
beam divided by the root-mean-square (rms) of the image resid- 
uals (i.e., after subtracting the CLEAN model). The peak flux 
densities of the images obtained from the self-calibrated visi- 
bilities range between 0.16 and 2.60 Jy at 8.4 GHz and between 
0. 16 and 1 .86 Jy at 15 GHz. The typical rms of these images are 
1 mJy beam" ' at 8 .4 GHz and 2 mJy beam" ' at 1 5 GHz. 

A total of 64 phase-referenced images were obtained at each 
frequency. These are two images for each pair of sources (i.e., 
image of source A phase-referenced to B and image of source B 
phase-referenced to A). We then discarded the images with dy- 
namic ranges lower than 10, because in these cases there is a high 
chance of confusion of the source with a spurious noise peak. 
Applying this cutoff to the dynamic ranges, 53 phase-referenced 
images were left at 8.4 GHz and 31 images at 15 GHz. The dy- 
namic ranges of the phase-referenced images are typically a fac- 
tor ~ 40 smaller than those obtained from the self-calibrated 
visibilities (see Figs. [T] and |2]). A first conclusion is that the loss 
of phase coherence strongly affects the dynamic range of the 
phase-referenced images, regardless of the calibrator-to-target 
separation. 

A similar conclusion was reported by Martf-Vidal et al. 
(2010) based on Monte Carlo simultations of phase-referenced 
observations under a turbulent atmosphere. These authors mod- 
elled the dynamic range of a phase-referenced image consider- 
ing the addition in quadrature of two sources of noise; one due 
to the thermal noise of the receiving system, cr,/,, and the other 
due to the atmosphere. o"a,. The latter was assumed to be equal 
to a given percentage of the source flux-density, i.e. a-gi - fm S . 
In this expression, S is the flux density and /„, is a factor that 
depends on the calibrator-to-target separation, 0, the observing 
frequency, v, and the on-source observing time Af. We notice 
that a similar expression is usually employed in the estimate 
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servations only depend on the observing frequency. The dynamic 
range is thus 



D = 



S 



(1) 
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The dynamic range of a phase-referenced image is Umited to 
a given value, Di, 



1 
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This limit is independent of both S and the sensitivity of the 
array. It is achieved when the flux density of the target source 
is much higher than the thermal noise of the interferometer The 
limiting dynamic range only depends on 6, v, and Af, and it can 
be several orders of magnitude smaller than the corresponding 
dynamic range due to the thermal noise of the receiving sys- 
tem, indicating that the atmosphere strongly limits the sensitiv- 
ity of the observations. Marti- Vidal et al. (2010) propose a phe- 
nomenological model for /«, based on their Monte Carlo simu- 
lations (see their Sect. 5.3). This model takes the form 



fa, = K- 



VAfsin (ff) 



(3) 



Fig. 1. Dynamic ranges (normalized to an observing time of 10 
hours) of the phase-referenced images as a function of distance 
to the calibrators. The error bars are proportional to the flux den- 
sities of the calibrators. Lines represent a model of the maximum 
achievable dynamic range (Eq.|2]i. 
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Fig. 2. Dynamic ranges (normalized to an observing time of 10 
hours) of the images obtained from self-calibrated visibilities. 
Stars are data at 8.4 GHz and diamonds at 15 GHz. This figure 
can be compared to Fig. [T] to see the effect of coherence losses 
in the phase-referenced images. 



of the (instrumental) systematic amplitude-calibration errors in 
VLA and VLBI observations. However, in our case the factor 
fat depends on the calibrator-to-target separation and on the to- 
tal on-source observing time, while the factors for the estimate 
of systematic amplitude-calibration errors in VLA and VLBI ob- 



where A' is a constant to be determined. 

We can compare this phenomenological model to our obser- 
vations by plotting the dynamic range of all our phase-referenced 
images as a function of calibrator-to-target separation. The im- 
ages obtained for each source pair correspond to data with 
slightly different on-source observing times. Therefore, we cor- 
rected (i.e., we normalized) the dynamic ranges of all the im- 
ages by applying the factor y/Nio/Ny, where A^,, is the number 
of visibilities for each image and A^io is the number of visibil- 
ities corresponding to a on-source observing time of 10 hours 
with the VLBA. We show in Fig. [1] the normalized dynamic 
ranges and the model resulting from Eq. |3] We find that for 
K ~ 12.4 hr"-'' GHz ' the model predicts the limiting dynamic 
ranges obtained with the VLBA at both frequencies, although 
the results at 8.4 GHz are of higher quality; the results at 15 GHz 
are noisier The dynamic ranges (also normalized to an observ- 
ing time of 10 hours) obtained from the self-calibrated images 
are shown in Fig.|2]for a comparison with those obtained from 
phase-referencing. The error bars in Fig. [1] are set to be propor- 
tional to the flux density of the calibrator. This way the reader 
has information in a single figure on the quality of the phase- 
referenced images and the quahty of the caHbrator visibilities. 

3.2. Peak flux-density 

We show in Fig.|3]the peak flux densities of the phase-referenced 
images relative to those obtained from the images correspond- 
ing to the self-calibrated visibilities as a function of distance 
to the calibrator. The systematics in the loss of flux density is 
clear for the 8.4 GHz data. The flux density recovered is about 
80% of the real flux density for separations of ~ 5 degrees, 
and slowly decreases to 50% for the largest separations. For 
the 15 GHz data, the scatter is larger and no robust conclusion 
can be obtained. However, a hint of saturation of the peak ratios 
around 40% can be appreciated at large source separations if the 
points corresponding to a separation of 14.84 degrees (for which 
peak ratios as high as 80% are obtained) are not considered 
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Fig. 3. Peak flux densities of the phase-referenced images nor- 
malized to the peak flux densities of the images obtained from 
the self-caHbrated visibilities. Lines represent the model given 
byEq.m 



(these points correspond to sources B 1803 H-784 and Bl 150H-812 
phase-referenced to each other). 

From the simulations reported in Martf-Vidal et al. (1201 Oi l, 
the peak flux-density of a phase-referenced image, Ppf,, relative 
to the real peak flux-density of the source, Psf, can be estimated 
as 



where ki is a constant to be determined. From the Monte Carlo 
simulations, the exponent k2 takes the value 1.87 ± 0.02. In Fig. 
[3] we also show the predictions of the model given by Eq. |4] 
for ki ~ 63. The model roughly predicts the behaviour of the 
8.4 GHz data, although for separations larger than ~ 15 degrees, 
the observed peak ratios clearly saturate at ~ 0.5, while the 
model predictions monotonically decrease. A similar conclusion 
is obtained for the 15 GHz data, where a saturation around 0.4 is 
appreciated. The saturation in the ratio of flux-density peaks at 
large calibrator-to-target separations is not modelled using Eq. 
|4] and could be due to the saturation of the power-spectrum of 
the tropospheric turbulences at large scales (see, e.g., Thomson, 
Moran, & Swenson ll991l l. which was not considered in the sim- 
ulations reported in Martf-Vidal et al. (1201 Oi l. This saturation of 
the power spectrum of the turbulences would stabilize the phase 
difference between target and calibrator for large separations and 
therefore enhance the signal coherence (and peak flux-density) 
of the target source. 



4. Conclusions 

We report how phase referencing affects the signal coherence 
(and the fidelity of the images) in VLBI observations at different 
frequencies (8.4 GHz and 15 GHz) and for different calibrator- 
to-target separations (from 1.5 to 20.5 degrees). We determined 
the loss of dynamic range and peak flux-density of the phase- 
referenced images and compared the results with the model pre- 
dictions given in Martf-Vidal et al. ( 1201 Oi l. 

The dynamic range of the phase-referenced images is 
strongly limited by the atmosphere at all frequencies and for 
all calibrator-to-target separations. The maximum achievable dy- 
namic range using the VLBA is given by Eqs. |2] and [3] with 
K ~ 12.4 (if V is given in GHz and At in hours). If the target 
source is not strong (as is usually the case), the thermal noise of 
the receiving system cannot be ignored and the dynamic range 
should be estimated with Eqs.[T]and|3] 

The flux-density (computed as the peak flux density of the 
phase-referenced image in units of the real peak flux-density 
of the source) decreases as the separation to the calibrator in- 
creases and is given by Eqs. |4]and[3] with A:i ~63andA;2 ~ 1.87. 
This model roughly predicts the peak ratios observed at 8.4 GHz 
for calibrator-to-target separations below ~ 15 degrees (the re- 
sults at 15 GHz are too noisy to reach a robust conclusion). For 
separations larger than 15 degrees, the observed peak ratios are 
higher than the model predictions, possibly due to a saturation 
in the power spectrum of the tropospheric turbulences at large 
scales, which was not considered in Martf-Vidal et al. (1201 Oi l. 

It is remarkable that for relatively small separations (below 
5 degrees), which are typical in many phase-referencing obser- 
vations, the flux-density loss can be as large as 20% at 8.4 GHz 
and 30 - 40% at 15 GHz (and even larger at higher frequencies, 
according to Eq. HJ. It must be also taken into account that the 
phase-referenced observations here reported were performed un- 
der good weather conditions and when the sources were close to 
their transits at nearly all stations (except Mauna Kea and St. 
Croix, see Martf-Vidal et al. 12008! for details on the observing 
schedule). Therefore even larger flux-density losses and lower 
dynamic ranges may be obtained when the observing conditions 
are far from optimal. 

However, we must also notice that the typical calibrator-to- 
calibrator cycle times in our observations were about 120 - 180 
seconds. For observations at 8.4 GHz, these cycle times are of 
the order of (and slightly shorter than) the coherence time of the 
signal (which takes typical values of 180 - 240 seconds in VLBI 
observations). Therefore it is not expected that a changing atmo- 
sphere may have introduced strong effects in the coherence of the 
phase-referenced visibilities at 8.4 GHz. At 15 GHz, the shorter 
coherence time (which takes typical values of 100 - 140 seconds 
in VLBI observations) might be short enough to imply an addi- 
tional phase degradation in the phase-referenced visibilities due 
to a changing atmosphere. Nevertheless, this issue is unlikely to 
significantly affect the results of our analysis, because the time 
evolution of the self-solutions of the antenna-based phase gains 
at 15 GHz are smooth, which indicates that the weather condi- 
tions were good enough to allow for a well-behaved phase con- 
nection between contiguous scans of each source. 
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